Dispersion spreading of biphotons in optical fibres and two-photon interference 
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We present the first observation of two-photon polarization interference structure in the second-order 
Glauber's correlation function of two-photon light generated via type-II spontaneous parametric down- 
conversion. In order to obtain this result, two-photon light is transmitted through an optical fibre and 
the coincidence distribution is analyzed by means of the START-STOP method. Beyond the experi- 
mental demonstration of an interesting effect in quantum optics, these results also have considerable 
relevance for quantum communications. 
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In many quantum communication, quantum computa- 
tion and quantum metrology experiments, entangled 
states of light are transmitted through optical fibres. For 
example, quantum key distribution through fibres has al- 
ready been demonstrated up to 100 km Because of 
these applications, a clear understanding of the effect of 
fibre propagation on the properties of entangled states 
is highly demanded. This understanding would be also 
very useful for the studies on the foundations of quantum 
mechanics 

In particular, it is known that an entangled two-photon 
state, a biphoton, spreads in time when propagating 
through a medium with group velocity dispersion. More 
specifically, the second-order intensity correlation func- 
tion of two-photon light gets broadened and at a suf- 
ficiently large distance, it takes the shape of the spec- 
trum Q , . It turns out that for biphotons propagating 
through a dispersive medium, the shape of the correlation 
function can manifest the effects of two-photon interfer- 
ence. 

In a typical experiment on observing two-photon in- 
terference for type-II spontaneous parametric down- 
conversion (SPDC) (Fig.l), pairs of orthogonally polar- 
ized signal and idler photons generated in the coUinear 
frequency-degenerate regime are sent to a non-polarizing 
beamsplitter (NPBS) followed, in both arms, by polar- 
izers (PI, P2) and photon counting detectors (Dl, D2). 
When one of the polarizers is rotated, the counting rate 
of the corresponding detector does not vary, but the coin- 
cidence counting rate varies with almost 100% visibility 
(polarisation fringes) 0. Another way to observe two- 
photon interference is to vary the phase between signal 
and idler photons (space-time fringes) ■ Such an exper- 
iment can also be interpreted as preparation of entangled 
two-photon Bell states in the two spatial modes after the 
beamsplitter 0- 

In addition to space-time and polarization fringes, two- 
photon interference can manifest itself in the modula- 
tion of the coincidence counting rate due to the displace- 
ment of one or both detectors in the near-field or far- 
field zone 01 or variation of the frequency registered by 



one or both detectors |lO|. However, to the best of our 
knowledge, there have been no experiments on observ- 
ing two-photon interference in the shape of the second- 
order intensity correlation function for SPDC, although 
a related experiment has been recently performed for a 
narrow-band two-photon source based on an OPO op- 
erating below threshold [l]| . The main reason is that 
normally, the width of the correlation function for SPDC 
radiation is on the order of hundreds of femtoseconds, 
and it cannot be resolved using the existing equipment. 
However, after a fibre of length about several hundred 
meters, the width of the correlation function grows to 
several nanoseconds, which enables the interference to 
be observed. 

In a setup shown in Fig.l, the state vector of the bipho- 
ton field at the output of the beamsplitter is 



l*) = 
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where F{Vl) is the spectral two-photon amplitude, Vl is 
the frequency shift, a^^ are photon creation operators 
in the horizontal and vertical polarization modes (de- 
noted hy a — H ,V ) and two spatial modes (denoted by 
i = 1,2); ujQ = i^p/'i,, Up being the pump frequency, and 
the terms corresponding to the case where both photons 
of the two-photon state go to the same port of the beam- 
splitter have been omitted. The difference between the 
group velocities of signal and idler photons in the non- 
linear crystal where SPDC occurs leads to the factors 
giiOro |-,y ^Yie two terms of Eq. (1), where tq = DL/2, 
with L being the crystal length and D = — ^ the 
inverse group velocity difference. To observe two-photon 
interference in type-II SPDC, this delay tq has to be com- 
pensated by means of a birefringent crystal. However, as 
one will see from what follows, compensation is not nec- 
essary in the presence of the fibre. 
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The coincidence counting rate can be calculated as 



Rc 



dt / dTG'-^\t,T), 



(2) 



where G^^^ {t, t) is the second-order Glauber's correlation 
function defined as 

G(2)(t,r) ^ {'f\Et\t)Ei~\t + T)Ei+\t + T)E^^\m), 

(3) 

e[^ (t) are positive- and negative- frequency field opera- 
tors at the photodetectors, integration in t is performed 
over the measurement time interval, and integration in 
T, over the coincidence circuit resolution window [l^. If 
SPDC is obtained from a cw pump (the stationary case) , 
G^^^ does not depend on t. 

One can show that two-photon interference man- 
ifests itself most explicitly in the difference between co- 
incidence counting rates for two cases. In the first case, 
both polarizers are oriented at 45 degrees to the horizon- 
tal direction. In the second case, polarizer 1 is oriented 
at 45° and polarizer 2, at —45° to the same direction. 

For these two configurations of the polarizers (45°, 45° 
and 45°, —45°), the correlation function takes the forms 



(r) = |F(r-ro)±F(T + ro) 



(4) 



where F{t) is the Fourier transform of the spectral am- 
plitude F{n) . 

Note that for type-II SPDC the amp litude F{t) has a 
rectangular shape with width 2tq 13*. Hence, if the de- 
lay To is not compensated, the two amplitudes in Eq. (3) 
do not overlap and the interference is absent: the coinci- 
dence counting rate, given by the integral of Eq. over 
T, is the same for (45°, 45°) and (45°, -45°) positions of 
the polarizers. 

If the biphoton beam, before entering the beamsplitter, 
passes through a sufficiently long fibre (with length z) 
|14| . the biphoton amplitude becomes [j] 



F(r) 



1 



V Airik" z 



(5) 



where fc', k" are, respectively, the first and the second 
derivatives of the fibre dispersion law k[ijj) at the fre- 
quency loq. Here, we assume that the group velocity dis- 
persion is the same for both polarizations. 

Since for type-H SPDC F(n) = sinc(Tol7) 0- 0, 



G^l>{e) ~ |e sinc(- 



To 



±e ^"0"/ sinc(- 



■ To. 



(6) 



where 9 = t — k'z is the shifted time and t/ = 2k" z/t^ 
is the tQ)ical width of the correlation function after the 
fibre lla . 



Taking into account that Tf » tq , one can rewrite 
Eq. (6) as 



^(2).^N sin^(g/T/) 

^(2)..^ Sin^(g/Tj)cOS^(0/T/) 



(7) 



Thus, in the presence of the fibre, two-photon interfer- 
ence can manifest itself in the shape of the second-order 
Glauber's correlation function. However, if the corre- 
lation function is integrated over all delays 9 (which is 
the case in most experiments on two-photon coincidence 
counting), the result will be the same for both config- 
urations of the polarizers. To see the structure given 
by Eq.s (jZj), one has to measure the correlation function 
with a sufficiently good time resolution. Two main rea- 
sons for the resolution reduction must be considered: the 
time jitter of the detectors (typically of the order of sev- 
eral hundreds of picoseconds), which sets the ultimate 
time resolution of the measurements, and the jitter con- 
tribution of the electronics (amplitude walk and noise) 
used for the coincidence detection. 

Coincidence detection by means of a coincidence logic 
gate circuit typically shows resolving time, for commer- 
cial available devices, of a few nanoseconds, which can 
be much greater than the spread r/ of the correlation 
function because of the propagation in a dispersive fibre. 

Another widely used coincidence detection technique, 
which will be further called the START-STOP method, 
involves direct measurement of the time delay between 
the photocount pulses of the two detectors by means 
of a time-to-amplitude converter (TAG), which converts 
linearly the time interval between the two input pulses 
(START and STOP) into an output pulse of a propor- 
tional amplitude. This analog pulse is forwarded to a 
multichannel analyzer (MCA) , which gives the histogram 
of the input pulse amplitudes corresponding to the prob- 
ability distribution for the time interval between the 
counts of the two detectors. One can show that in 
the limit of small photon fiuxes, this time interval distri- 
bution coincides in shape with the second-order intensity 
correlation function. The resolution of such technique 
could be on the order of one picosecond if the detector 
jitter time were negligible. The expected experimental 
coincidence distribution, in the case of negligible jitter 
time contribution from the experimental apparatus, was 
calculated numerically (Fig. 2a) for the following param- 
eters: k" ^ 3.2 • 10-28sVcm, fibre length z = 250m, 
D — 1.5ps/cm, crystal length L = 0.05cm. The cor- 
responding FWHM of the correlation function after the 
fibre is 1.2 ns (r/ = 0.43ns). 

One can see that for the polarizers set at (45°,— 45°), 
destructive interference occurs at the center of the peak. 
On the other hand, for the polarizers set at (45°, 45°), 
there is constructive interference at the centre of the 
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peak. High- visibility polarization interference can be ob- 
served in this case by selecting only several channels of 
the MCA corresponding to the center of the peak and 
rotating one of the polarizers. 

If we take into account the time jitter of the exper- 
imental setup, the interference structure gets smeared. 
In Fig. 2b we show the numerical evaluation of the coin- 
cidence distribution with an account for the time jitter 
of the experimental apparatus, which is the convolution 
of the peak plotted in Fig. 2a with a 750 ps (a typical 
jitter for two commercial APD detectors) FWHM Gaus- 
sian function. One can see that in the case of Fig. 2b, the 
change in the number of coincidences observed at zero 
delay for (45°,45°) and (45°, -45°) orientations of the 
polarisers corresponds to only 45% visibility, unlike in 
the case of Fig. 2a, where almost 100% visibility is ob- 
served. 

It is worth hinting to what happens if the e-o delay tq 
is compensated partly, by means of a birefringent mate- 
rial introducing an opposite-sign delay. Here [l^ we only 
mention that the period of the sine and cosine functions 
in Eqs. Q increases above the width of the envelope and 
thus the integral numbers of coincidences corresponding 
to the two polarizer configurations differ. As a result, 
there is nonzero visibility even in the absence of the cor- 
relation function time selection. For complete compen- 
sation, the visibihty would be 100%. 

In order to experimentally observe two-photon inter- 
ference in the shape of the correlation function of bipho- 
tons, we have generated biphoton pairs via SPDC by 
pumping a type-II 0.5 mm BBO crystal with a 0.1 Watt 
CW Ar"*" laser beam at the wavelength 351 nm in the 
coUinear frequency-degenerate regime. No compensating 
birefringent material was used after the crystal. Then 
the pump beam was eliminated by a high-reflectivity UV 
mirror plus an anti-UV cutoff filter and the SPDC ra- 
diation was coupled to a 250 m - long single-mode fi- 
bre with FMD= 4/Lt by a 20x microscope objective lens. 
In order to increase the efficiency of coupling SPDC ra- 
diation into the fibre, the pump was focused into the 
crystal using a lens with the focal length 13.5 cm. Be- 
cause of the rapid polarization drift in the fibre, a fi- 
bre polarization controller was introduced at the output. 
After the fibre, the biphoton pairs were addressed to a 
nonpolarizing beamsplitter and two photodetection ap- 
paratuses (consisting of polarisers, red-glass filters and 
APDs). The photocount pulses of the two detectors, af- 
ter passing through delay lines, were analyzed by means 
of the START-STOP method, and the second-order cor- 
relation function, G^'^\9)^ was observed at the MCA out- 
put. The FWHM of the coincidence peak in the absence 
of the fibre was measured to be 0.75 ns, a value substan- 
tially determined by the APD time jitter. The resolution 
determined by the MCA channel width, in our case 2.5 
ps, is much smaller and therefore negligible. 

If the biphotons propagate through the fibre, the 



FWHM grows up to 1.2 ns. The peak observed without 
polarization selection is shown in Fig. 3a. In the pres- 
ence of the polarizers, the shape of the peak changes 
(Fig. 3b, c). For the (45°;— 45°) settings of the polarizers, 
destructive interference in the middle of the peak leads 
to a shape in agreement with the calculated one (Fig. 2b). 
For the polarizers set at (45°;45°), as expected, there is 
a maximum at the center of the peak. If only the central 
part of the peak is selected, the observed variation in the 
number of coincidences due to the change in the polar- 
izer settings corresponds to 35% visibility of polarization 
interference. 

In summary, we have observed effects of two-photon in- 
terference in the shape of the second-order Glauber's cor- 
relation function for biphoton light generated via type-II 
SPDC and transmitted through an optical fibre. It is 
worth noticing that our results show how, for a suffi- 
ciently long fibre (hundreds of meters), a Bell state can 
be generated without compensating for the e-o delay be- 
tween the signal and the idler photons. In this partic- 
ular method of Bell states preparation, distinguishabil- 
ity of the interfering two-photon amplitudes is erased 
by spreading them in time and by selecting coincidences 
within a suitable time window. Even if this method could 
appear more complicated than the traditional one, based 
on using additional birefringent elements, it is interest- 
ing as a conceptually different way for restoring entangle- 
ment. In this sense it must be noticed that selection of 
the coincidence time window is similar to the frequency 
selection, which can be achieved by using narrow-band 
filters in some experiments on Bell-state preparation 
since the group- velocity dispersion of the fibre ac- 
tually performs the Fourier transformation of the bipho- 
ton spectrum into the correlation function How- 
ever, this does not mean that the observed interference 
is present in the spectrum of the down-converted pho- 
tons, since the observed shape of the correlation function 
substantially depends on the polarizer settings in both 
arms, which, in the general case, are different. Another 
distinction 'l.^ from the experiments where narrow-band 
filters are used for the Bell states preparation consists of 
the fact that, in our experiment, one can achieve high vis- 
ibility by selecting not necessarily the central part of the 
two-photon peak but any symmetric side parts where the 
coincidence counting rates differ for different polarizer 
settings (see Fig. 2). Furthermore, one should notice that 
these effects inevitably occur in any quantum commu- 
nication experiment where biphoton light is transmitted 
through optical fibres, and if the fibres are sufficiently 
long and temporal selection of the coincidences is pro- 
vided, there is no need for o-e delay compensation. 

Finally, we would like to mention that the measure- 
ment of the observed broadening of the second-order 
Glauber correlation function could eventually be used for 
the evaluation of th e op tical fibre chromatic dispersion 
(in analogy to Ref. flTf). 
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Optical fibre 




FIG. 1: A setup for observing two-photon polarization inter- 
ference for type-II SPDC. The fibre is introduced for spread- 
ing the correlation function and observing its interference 
structure. 





FIG. 2: Time delay distribution of coincidences calculated for 
k" = 3.2-10~^**sVcm, z = 250m, D = 1.5ps/cm, L = 0.05cm, 
and the polarizers oriented at (45°,— 45°) (solid line) and at 
(45°, 45°) (dashed line) (a) for an ideal measurement setup 
with negligible time jitter; (b) for a measurement setup with 
a time jitter of 750 ps. As the mutual orientation of the 
polarizers changes from parallel to orthogonal, the change in 
the coincidence number at zero delay corresponds to 100% 
visibility in the (a) case and to 45% visibility in the (b) case. 
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FIG. 3: Observed time delay distribution of coinci- 
dences: (a) without the polarizers; (b) with polarizers set at 
(45°, 45°), which corresponds to measuring G^l\9); (c) with 
polarizers set at (45°, —45°), which corresponds to measuring 
G^_^\o). Destructive interference at the center of the coinci- 
dence peak is clearly seen. As the parallel polarizer configura- 
tion (b) is switched to the orthogonal polarizer configuration 
(c) , the change in the coincidence number at the center of the 
peak corresponds to approximately 35% visibility. 



